Evidence for Static Magnetism in the Vortex Cores of Ortho-II YBa2Cu3O6.50 
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Evidence for static alternating magnetic fields in the vortex 
cores of underdoped YBa2Cu306+a; is reported. Muon spin 
rotation measurements of the internal magnetic field distribu- 
tion of the vortex state of YBa2Cu306.5o in applied fields of 
H = l T and H = 4T reveal a feature in the high-field tail of 
the field distribution which is not present in optimally doped 
YBa2Cu3 06.95 and which fits well to a model with static mag- 
netic fields in the vortex cores. The magnitude of the fields 
is estimated to be 18(2) G and decreases above T = 10 K. 
We discuss possible origins of the additional vortex core mag- 
netism within the context of existing theories. 



The application of a large magnetic field to a supercon- 
ductor drives part of the sample, the vortex cores, into a 
"normal" state. While the physics of these vortex cores 
in conventional superconductors is generally thought of 
as metallic, the vortex cores of high temperature super- 
conductors may offer insight into more unusual low tem- 
perature properties. A new class of theories has predicted 
that magnetism may be induced near/inside the cores of 
vortices by the application of a magnetic field j^]- |Q. For 
example in Zhang's SO{5) theory, static antiferromag- 
netism (AF) should appear in regions where the super- 
conducting (SC) order parameter is suppressed, and as 
a consequence, vortices in underdoped superconductors 
should be magnetic Similarly, Lee and Wen have 

predicted a staggered flux phase (SFP) of orbital currents 
in the vortex cores of underdoped cuprates Below 
Tc, the system is locally a d-wave superconductor away 
from the vortex cores, but has a SFP "frozen" inside the 
vortex cores. Consequently, they predict the appearance 
of quasi-static alternating magnetic fields of order 10 G 
in the vortex cores. Finally, Zhu and Ting have devel- 
oped the Hubbard model with an on-site repulsion term 
and shown that AF-like spin density waves (SDW) can 
appear in the vortex cores 1^, 

The possible coexistence of magnetism and supercon- 
ductivity in the cuprates has been a key issue since 
the discovery of high- Tc superconductors. In zero 
field, microscopic coexistence of superconductivity, mag- 
netism, and spin-glass behavior has been detected with 
/iSR in the underdoped region near the boundary be- 



tween SC and AF [p|, in Ca-doped YBa2Cu306-i-a; [§, 
La2_5Sr5Cu04 ||9|,|10[, and underdoped and optimally 
doped YBa2Cu306+3; pl] , p2| . There is also recent exper- 
imental evidence from neutron scattering for anomalous 
magnetism in underdoped YBa2Cu306+x |ll3| , |l4| in zero 
applied field. Field-induced low- frequency magnetic fluc- 
tuations in La2-5Sr5Cu04 have been reported by Lake 
|l| and in YBa2Cu306+:r by Mitrovic |l|l. Vaknin |l|] 
found a possible signature of AF cores in optimally-doped 
YBa2Cu306-i-a;, while Katano report sharp incommensu- 
rate neutron scattering peaks in La2-5Sr5Cu04 that are 
enhanced in a magnetic field [ p^ . 

In this paper, we report a search for static magnetism 
in the region of the vortex cores of underdoped Ortho- 
II YBa2Cu306.5o- Signiflcant improvements of the fits 
to the muon spin precession signal in the vortex state 
are obtained using a model of the vortex lattice with an 
additional alternating magnetic field of 18(2) G, whose 
magnitude decreases away from a vortex core center on 
the length scale of the coherence length. As a control, 
similar fits were made on data from the conventional su- 
perconductor NbSe2 in which case no improvement in the 
fits was observed, as expected. 

Muons are an excellent probe of magnetism in super- 
conductors. As described elsewhere ||l^, the implanted 
spin polarized muons stop randomly on the length scale 
of the vortex lattice and precess at a rate proportional to 
the local magnetic field, thus providing a direct measure 
of the local field distribution n(B). The observed asym- 
metric magnetic field distribution in a superconductor is 
characteristic of a lattice of magnetic vortices in a type- 
II superconductor. In Ginzburg-Landau (GL) theory, a 
vortex core's size is determined by the applied magnetic 
field H and the in-plane GL coherence length S^ab, while 
the magnetic field decays away from the vortex core over 
a length scale given by the GL ab in-plane penetration 
depth Xab- The magnetic field distribution n{B) can be 
calculated from the spatial distribution of the magnetic 
field H: 

where = 2^1^G'^{1 + b*)[l - 2b{l - b)% Ki{u) is a 
modified Bessel function, G is a reciprocal lattice vector 
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of the vortex lattice, b = H/Hc^ is the reduced field, $o 
is the flux quantum, and S is the area of the reduced unit 
cell for a hexagonal lattice. 

A phenomenological model of n{B) that includes alter- 
nating fields in the vortex cores can be created by adding 
a term to Eq. |lj 

B'{r) = B{r) + (-l)(^+s')/'^Me-(''/«"^)', (2) 

where M is the amplitude of the alternating fields in the 
vortex core that alternates in sign between neighboring 
crystal lattice sites (a is the distance between sites) . Our 
model assumes that the amplitude of the fields in the vor- 
tex cores decays away from the vortex core center over 
the same length scale S^ab as the SC order parameter in- 
creases. We note: 1) that Eq. ^strictly holds only in the 
unit cell of the vortex lattice, in the limit ^ab/Lo 1, 
where Lq is the distance between vortex cores. This en- 
sures there is negligible overlap between alternating fields 
of neighboring vortices. 2) Because the vortex and crystal 
lattices are incommensurate, the parameter a has been 
varied between the Cu-Cu distance in the planes and a 
much smaller value that allows the muon to experience all 
alternating fields from ~M to M. The results presented 
here are insensitive to such a change. 

Figure |l| shows the effect of the vortex core alternat- 
ing fields on n{B). The dashed line shows n{B) without 
alternating fields (M = G) , for typical values A = 2450 
A and C = 35 A in an applied field of = 4 T. The 
high-field cutoff corresponding to the magnetic field at 
the vortex core center is visible at b. The modified n{B) 
due to the introduction of the fields of magnitude M — 15 
G (solid line) shows a double step in the high field tail 
at a and c. The lower step (at a) is caused by additional 
fields in the vortex core antiparallel to the applied mag- 
netic field H, whereas the upper step (at c) originates 
from sites where the additional fields are aligned with 
the applied field. The spatial variation of the magnetic 
field near the vortex cores is shown in the inset. The 
line corresponds to B{r) from Eq. |l| without alternating 
fields, while the squares show B'{r) from Eq. |^ with the 
fields. The van Hove singularity in n{B) at a is caused 
by the nearly flat field distribution in the inset along the 
lower line of square points. These points correspond to 
fields ahgned opposite to the appHed field. 

With a model for n{B), a theoretical function of the 
muon spin precession signal P{t) =Px+iPy is constructed 
from Eq. |^ or || by calculating the local magnetic field at 
points on a lattice in the reduced unit cell and summing 
the contributions from all sampling points: 

/"OO 

Pa:{t) ^ P{Q)G{t) n{Lu)cos{ujt + 5)duj, (3) 
Jo 

where uj = 'y^.B, G{t) is a gaussian relaxation function 
that accounts for additional broadening due to nuclear 



dipole moments and lattice disorder pi| , and 27r7^ = 
135.54 MHz/T. 

In our fits to the muon spin precession signal, n{B) 
parameters £,ab, ^ab, M, and H were allowed to vary. It 
has been shown previously that Eqs. |l| and [s] can be 
used to model the muon spin precession signal in NbSe2 
p2| , p^ and YBa2Cu3 06+a; [ p6|j25| superconductors. In 
this work, constraints on the fits of Eq. ^ were applied 
to the range of allowed values of £_ab- We find that at 
the highest applied fields (if = 4 T), ^ab « 25 A, a 
value smaller than expected from previous /LtSR measure- 
ments of YBa2Cu3O6.60- Therefore, in the fits, £,ab is 
constrained to ^ab> 35 A at H = AT. Unconstrained fits 
yield larger values for the alternating field magnitude. 

The /xSR experiments were performed on the M15 sur- 
face muon channel at TRIUMF using a horizontal gas- 
cooled cryostat and a high-field (7.5 T) high timing res- 
olution spectrometer (BELLE). All data were taken un- 
der conditions of field cooling, and high statistics runs 
were made with approximately 30 million muon decay 
events for each value oi B,T. A mosaic of high pu- 
rity YBa2Cu306.5o single crystals were used. They were 
grown in BaZrOa crucibles, and detwinned below 200 C. 
Ortho-II ordering was achieved by low temperature an- 
nealing for 7 days |^ . The oxygen content of the sample 
was chosen due to its stoichiometry and high degree of 
crystalline order, proximity to the AF/SC boundary, and 
the oxygen ordering of the Ortho-II phase. 

/iSR measurements performed in zero magnetic field 
in YBa2Cu306.5o show no internal magnetic fields. Fur- 
thermore, the magnetic field dependence of the skewness 
and second moment of the field distribution shows no 
signature of a transition to a 2D pancake state ||2^ at 
r = 2.2 K, where n{B) is no longer described by Eqs. |l| 
and ^. In a field of 4 T, an asymmetric field distribution, 
the signature of a well-ordered lattice, is visible up to 
r = 30 K. Further details of the vortex phase diagram in 
YBa2Cu306.5o will be published elsewhere. 

Figure ^ shows a fit to the muon spin precession signal 
in YBa2Cu306.5o assuming the phenomenological field 
profile of Eq. ||. Notice the excellent fit of the theoretical 
line to the data, and the small uncertainty in the data 
due to the high statistics. The significant improvement 
in the fits of our model with alternating magnetic fields is 
shown in Figure ^ Fig. ^(a) shows the deep minimum in 

/NDF (number of degrees of freedom) achieved by the 
addition of fields of magnitude A/ = 15 G. Figs, ^(b) and 
(c) show that a similar effect is not found in optimally 
doped YBa2Cu3 06.95 (at H — 6 T) or in the conventional 
superconductor NbSe2 (at H = 0.7 T), which is generally 
well described by BCS theory. The absence of this effect 
in these materials indicates that the improvement of the 
fits is not simply due to the addition of another parame- 
ter and that the feature in YBa2Cu306.5o is not present 
in optimally doped YBa2Cu3 6.95 or conventional super- 
conductors. The curves in Fig. @ are guides to the eye. 
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The improvements in the fit can be seen more clearly 
in the frequency domain. FFT's of data from Fig. || and 
the corresponding fits are shown in Fig. ^. The triangles 
are an FFT of the data, whereas the full (dashed) line is 
a FFT of the fit with (without) vortex core fields. Fig. 
^(b) shows the high field tail of the Fourier transforms 
in greater detail. The improved fit is clearly visible in 
the high field tail, despite the broadening introduced by 
Fourier transforming finite data. Additional weight is 
present at fields higher than the high field cutoff found 
using Eq. 0. Similar improvements in the fits were found 
at other temperatures in applied fields oi H — 4 T and 
H=l T. The amplitude of the magnetic fields M in the 
vortex core center decreases gradually with temperature 
above T= 10 K. Beyond T=20 K, no clear signature of 
the fields could be seen. 

In our experiment, the muon is sensitive to the c-axis 
component of the local magnetic field at the muon site. 
Estimating the size and direction of the effective mag- 
netic moment responsible for these alternating fields from 
these results is not easy. However, Hsu et al. calculated 
that the geometric factor due to the displacement be- 
tween the muon site and the center of a CuO plaquette 
would be of order 1 [^, assuming the muon bonds to 
the oxygen above the CuO planes. Hence measurement 
of 18 G static fields at the muon site in the center of the 
vortex core agrees quantitatively with theoretical predic- 
tions Various other forms of the alternating field 
term in Eq. |^ have also been tried. Replacing the gaus- 
sian function in the alternating field term by an expo- 
nential yields a worse fit (although is still a minimum 
at non-zero Af), while fitting with a second length scale 
for the magnetic field amplitude decay function (instead 
of ^afc) provides only a marginal improvement in the fit. 

We note that there may be a number of other possible 
interpretations of the anomalous n{B) we have measured. 
For example, little is known about vortices in underdoped 
high-Tc superconductors. The d-wave nature of the vor- 
tices or changes in the pairing symmetry may affect n{B). 
However, to our knowledge no theories currently exist 
that predict a change in the high-field tail due to such 
effects. Another possibility is that the vortex lattice is 
characterized by two domains that coexist in the sam- 
ple, one a well-ordered vortex flux lattice and the other 
a disordered flux lattice. The high quality of the sample 
makes the likelihood of such coexistence small Thus, 
we are led to the conclusion that the most likely interpre- 
tation is the presence of alternating fields in the vortex 
cores. As noted above, numerous theories predict their 
existence, and experimental evidence of fluctuating AF 
fields in the vortex cores of La2-5Sr5Cu204 has been re- 
cently reported in neutron scattering measurements. Our 
main result is precisely the signature one would expect in 
the field distribution of a superconductor with additional 
alternating magnetic fields in the vortex cores. 

In conclusion, we have shown that a phenomenological 



model of the field distribution with alternating magnetic 
fields in the vortex cores fits the muon spin precession sig- 
nal in YBa2Cu306.5o in high magnetic fields significantly 
better than a model without such vortex core fields. We 
find the amplitude of these fields to be 18(2)G at the 
muon site. The origin of the anomalous magnetism re- 
ported here may be consistent with Zhang's 5*0(5) spin 
magnetism or various t — J Hubbard models of orbital 
magnetism. Further measurements at intermediate dop- 
ing between x — Q.h and 0.95 are in progress. This mea- 
surement is of significant importance to a class of theories 
of high-Tc superconductors that predict spin or orbital 
magnetism in underdoped cuprates. 



ACKNOWLEDGMENTS 

We wish to gratefully acknowledge discussions with Ian 
Affleck and the help of Kim Chow, Mel Good, Bassam 
Hitti, Don Arseneau and Syd Kreitzman and of the tech- 
nical staff at TRIUMF. Correspondence should be ad- 
dressed to R.I. Miller (miller@physics.ubc.ca). 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 

[s: 

[9 

[lo; 
[11 

[12 
[13 

[14 
[15 
[16 
[17 

[18 
[19 

[20 
[21 
[22' 
[23^ 

[24' 
[25^ 



S.C. Zhang, Science 275, 1089 (1997). 

D. Arovas et al, Phys. Rev. Lett. 79, 2871 (1997). 

X. Wen and P. A. Lee, Phys. Rev. Lett. 7 6, 503 (1996). 
P.A. Lee and X. Wen, |cond-mat/000841S|. 
J.X. Zhu and C.S. Ting, |cond-mat/0105417i 
M. Franz and Z. Tesanovic, |cond-mat/0002137 , 

E. Demler, S. Sachdev, Y. Zhang, Phys. Rev. Lett. 87, 
067202 (2001). 

J.H. Brewer et al, Physica C 162, 157 (1989). 

Ch. Niedermayer et al, Phys. Rev. Lett. 80, 3843 (1998). 

S. Wakimoto et al, Phys. Rev. B 63, 172501 (2001). 

J. Sonier et al, Science 292, 1692 (2001). 

J. Sonier et al, zond-mat/010847c . 

H.A. Mook, Pengcheng Dai, F. Dogan, 



cond-mat/0102047 



Y. Sidis et al, |cond-mat/0101095 
B. Lake et al. Science 291, 1759 (2001). 
V.F. Mitrovic et al, Nature 413, 501 (2001). 
D. Vaknin, J.L. Zarestky, L.L. Miller, Physica C 329, 
109 (2000). 
S. Katano et al 
J.E. Sonier, J.H, 
724, 769 (2000). 

A. Yaouanc et al, Phys. Rev. B 55, 11107 (1997). 
Here, G{t) = exp{-{atf /2). See Ref. (l9|. 
J.E. Sonier et al, Phys. Rev. Lett. 79, 1742 (1997) 
See discussion on vortex state phases in Ref. [llSl 
G. Blatter et al, Rev. Mod. Phys. 66, 1125 (1994) 
R.I. Miller et al, Phys. Rev. Lett. 85, 1540 (2000). 
J.E. Sonier et al, Phys. Rev. Lett. 79, 2875 (1997) 
Sonier et al, Phys. Rev. B 55, 11789 (1997). 



Phys. Rev. B 62, R14677 (2000). 
Brewer and R.F. Kiefl, Rev. Mod. Phys. 



and 



J.E. 



3 



[26] J.E. Sonier et al, Phys. Rev. Lett. 72, 744 (1994). 
[27] R. Liang, D.A. Bonn, and W.N. Hardy, Physica C 336, 
57 (2000). 

[28] T.C. Hsu, J.B. Marston, and L Affleck, Phys. Rev. B 43 
2866 (1990). 




0) 

E 
E 

< 




FIG. 2. Muon Polarization Signal Pj:{t) in YBa2Cu3O6.50 
in an applied field of H — 4 T and at T = 5 K displayed in a 
rotating reference frame. The line is a fit to Pa;(t) using Eqs. 
^ and |[. 
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FIG. 1. Magnetic field distribution n{B) in superconduc- 
tors in an applied field of iif = 4 T, with and without vortex 
core alternating fields. Solid line: n{B) calculated from Eq. |^ 
with alternating fields in the vortex cores (M = 15 G); dashed 
line: field distribution without alternating fields calculated 
from Eq. |l| with identical values for ^ and A. Notice that sig- 
nificant changes occur only in the high field tail. See text for 
further details. Inset: B{r) along a line between two vortex 
cores. The solid line shows B{r) without alternating fields, 
whereas the squares show B{r) with M = 15 G alternating 
fields in the vortex core. The magnitude of the alternating 
fields is assumed to decay away from the vortex core center 
on the length scale of ^ab- 



4 



3.0 



(a)YBaCuO,„ 

^ ' 2 3 6.50 

2.5- / 




1.0- 

10 20 30 
AF Field Amplitude, M, (G) 

FIG. 3. The best per degree of freedom versus the 
amphtude of the alternating field, M: (a) YBa2Cu3O6.50 at 
//■ = 4 T and at r = 5 K (solid squares); (b) YBa2Cu3 06.95 at 
H = 6 T (solid nablas); and (c) for NbSe2 at T (open 

squares). The solid curves are guides to the eye. 
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FIG. 4. (a) Fast Fourier transform of /iSR data (triangles) 
and fit (solid line) in Fig. |^ and of the best fit to Eqs. ^and 
^without alternating fields (dashed line), (b) The high-field 
region in greater detail. Notice the improved fit in the high 
field tail due to the alternating fields, corresponding to the 
vortex core area. 
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